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I. INTRODUCTION 



Gravitational collapse with formation of black holes is one of the main problems of general 
relativity. Choptuik |L] discovered a critical behavior of the gravitational collapse of massless 
scalar field by a numerical study. His result can be summarized as follows. Let {Si(p)} be 
one-parameter families of solutions of Einstein's equations for spherically symmetric massless 
scalar field minimally coupled to gravitational field, where p is a parameter which smoothly 
specifies the initial value of the scalar and gravitational fields and i is an index which specifies 
one of the several one-parameter families. For each family Si(p) there is a critical value p* 
of the parameter p such that solutions Si(p > p*) contain a black hole while Si(p < p*) 
do not, which are referred to as supercritical and subcritical, respectively. He found that 
near-critical solutions satisfy the following: (1) the strong field region is universal in the 
sense that it approaches the identical spacetime for all families (i.e. for all i), (2) the strong 
field region has a discrete self similarity, and (3) for supercritical solutions the black hole 
mass is proportional to a certain power of the deviation of the parameter p from the critical 
value p* as 



and the critical exponent (3 is universal, i.e. independent of i. The numerical value of the 
(3 was given as (3 ~ 0.37, and was speculated to be 1/e. Abrahams and Evans @] found 
the similar phenomena in axisymmetric collapse of gravitational wave with almost the same 
value of the critical exponent (3 ~ 0.38. Evans and Coleman || also found the similar 
phenomena with (3 ~ 0.36 in spherically symmetric collapse of radiation fluid, in which case 
the self similarity is not discrete but local one. The difference seems to come from that in the 
first two cases the matter propagates as waves while in the last case it does not. There might 
be some typical wavelength in the first two cases. It should be noted that these calculations 
may suggest that the meaning of the universality in (3) can be extended to independence 
of the critical exponent from details of systems, though in the first suggestion it meant its 
independence of initial data. 
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The phenomena described above are very striking, but there are almost no theoretical 
account for the essense of them, especially for the value of the critical exponent (3 and 
the mechanism to accquire the self similarity and the universality. Oshiro, Nakamura and 
Tomimatsu Husain, Martinez, and Nunez [||], and Brady [0 studied the critical behavior 
of scalar field collapse using the exact self-similar solution. They gave the the critical 
exponent (3 < 0.5. Their argument on (3 was not very quantitative, which was mostly due 
to the difficulty of estimating the form of the metric in the asymptotic region. 

If the critical behavior has very universal feature we may expect the similar behavior in 
seemingly very different models, which may enable us to investigate it in a simple model. 
On the other hand, the model should have some essential properties which the above models 
have in common. First, the spacetime should be asymptotically flat. Second, there must be 
energy transportation to the null infinity. This is because otherwise the spacetime necessarily 
ends up with a black hole with the mass of the initially prepared matter whenever a black 
hole forms (apart from the possibility that a star forms outside a black hole of small mass) 
so that there would be no criticality, which is the case of the Lemaitre-Tolman-Bondi 
spacetimes. Third, the matter content was "massless" in all of the above numerical studies. 
It is preferable that our model should also have the similar property. 

A simple model possessing these properties is the collapse of a thin shell coupled with 
outgoing null flux. We will find such a model which exhibits the critical behavior in black 
hole formation. We treat it analytically, which may help find the mechanism of the critical 
behavior. We shall find the critical exponent and its dependence on the details of the 
systems. 

Sec. H is a brief review of the treatment of a thin shell in general relativity. In sec. [I 1 1 



we present the model of thin shell coupled with null flux. We find its critical behavior in 
sec. [TV[ In sec. [V], approximate solutions of the dynamics of the thin shell is derived and 
the value of the critical exponent is calculated by using them. The exactness of the value 
is rigorously demonstrated in Sec. [VJ. Sec. |V11| is devoted to summary, and Sec. |V1U| to 
discussion. 



II. GENERAL TREATMENT OF DYNAMICS OF A THIN SHELL 



We briefly summarize the prescription of dynamics of an infinitely thin shell using metric 
junction which was first formulated by Israel ]7|. Let S be a compact 2-surface which 
represents the thin shell. The shell's trajectory is given by two embedding of S x R into the 
exterior and the interior spacetimes V + and V-, which are supposed to be the solution of 
Einstein's equation. Let n a denote the outward-pointing unit normal vector to the trajactory 
and h ab denote the induced metric on the trajectory. The nontrivial components of Eintein's 
equation on the shell are 

G ah n a n b = 8nT ab n a n b , (2.1) 
G ab n a h b c = 87rT ab n a h b c , (2.2) 

which are equivalent to 

{3) i? + K ab K ab -K 2 = -16nT ab n a n b , (2.3) 
D b K\ - D C K = 8nT ab n a h b CJ (2.4) 

where ( 3 'R denotes the 3-curvature of the trajectory and D denotes the covariant derivative 
with respect to the metric h ab . The difference of the above equations between V+ and V_ is 
interpreted as due to the energy-momentum of the shell: 

K ab S ab =[T ab n a n% (2.5) 
D b S b c = -[T ab n a h b c ], (2.6) 

where K ab = (1/2)^++^), [A] = A+-A-, and S ab is defined as S ab = [l/8n)[h ab K-K ab ] 
which is interpreted as the energy-momentum of the shell. 

We assume that the shell has the surface energy density a and tension (, i.e. , 

Sab = ou a u b - C{h ab + u a u b ), (2.7) 

where u a is a timelike unit vector along the shell's trajectory. We also assume that the shell 
is spherically symmetric. Then the parts tangent to and normal to the vector u a of eq. ( |2.G| ) 
are 
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u b D b a + ^L^(a -() = [T ab n a u% (2.8) 
ri 

u h D b u c = 0, (2.9) 

respectively, where R is the area radius of the shell. Eq. fl2.8|) is the equation for the energy 
flow. Eq. ( |2.9|) is the equation of motion of the fluid particle of the shell, which indicates 
that each particle moves along the geodesic on the three-dimensional hypersurface of the 
trajectory of the shell. 

From the definition of S ab one has 

<j = --L{u c u d + h cd )[V c n d ] 
1 \n d V d r] 

= - s L_±i (2 . 10) 

where r is the area radius of the symmetric sphere and [n^V^r] is evaluated at r — R. In 
terms of the shell's proper mass M. = 4iraR 2 it can be written as 

- 4r = K v ^]- ( 2 - n ) 

ri 

The equation for the tension ( can be also derived, but it is not an independent equation 
in the model presented below. The model also satisfies eq. (|2.5| ). Eqs. (|2.8| ) and ([2.1 0|) 
determine the motion of the shell. 



III. COLLAPSE OF A THIN SHELL COUPLED WITH NULL FLUX 

We assume that the outer side of the shell is coupled with outgoing null flux which may 
transport the energy of the shell to infinity and the inner side of the shell does not have 
any interaction. The exterior spacetime is the Vaidya spacetime. The interior spacetime is 
assumed to be flat, i.e. , is the Minkowski spacetime. We can write the line element as 

ds 2 = -(1 - 2$)dn 2 - 2dudr + r 2 dQ 2 , (3.1) 

where u is the retarded time, r the area radius, and dQ 2 the line element on a unit sphere, 
and 
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M(u)/r (exterior) ^ 



A similar situation is considered in || in the context of the creation of man-made universes. 
Note that M(u) gives the Bondi-Sachs mass of the null infinity at u. 

Let t denote the proper time of the shell. Let r = R(t) represent the trajectory of 
the shell. Note that R is the shell's radius. The r-component of the vector u a is written 
as u r = R, where a dot denotes the differentiation with respect to r. It follows from the 



normalization condition that the w-component is written as u u = (y R 2 + 1 — 2$ + R) 1 



(Hereafter $ denotes $+(«, R)-) Then the vector n is given as n a = s(—u u , \J R? + 1 — 2$), 
where s is plus/minus unity if r is increasing/ decreasing in the direction of n a . Eq. (|2.10|) 
yields 

R 2 + 1 = f, (3.3) 

where 

,l(M) + *(»Y l . (3.4) 



2 V R) \R 

The outgoing flux on the shell is 



, sMu u /A-nR 2 (exterior) 
[T ab n a u b ] = { (3.5) 

(interior) 



where we used M = M, a u u . Then eq. ( p.8|) yields 

2R. M 1 . . 

a + —(a-C) = s = r . (3.6) 

R 47rR ^R 2 + l-2$ + R 

We need two more conditions to specify the system completely. One is the equation of 
state of the shell and the other is the coupling of the shell and the spacetime. We assume 
that the former is given by 

C = ~W (3-7) 
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with rj a constant. For simplicity we assume that the latter is given by 



(aR^y = 0, (3.8) 



aR 1 -" = C747T, (3.9) 

where u and C are constants. In the u = — 1 case the mass M of the shell conserves, while 
uj = 1 is said to hold in the case of false vacuum bubbles ||. So our assumption is a slight 



generalization of them in the parameter space. From eqs. ( |3.4| ) and (|3.9|) one has 

f = ET u - 1 M + R a /2, (3.10) 

and eq. (|3~E|) yields 



M = spRR UJ (\/R 2 + l-2<S> + R), (3.11) 



where p = 2r} + u + l, M = CM,R = CR, and $ is given in eq. Q. Eqs. (gl| and 
are the basic equations in the subsequent discussions. 

IV. MODEL WITH CRITICAL BEHAVIOR 



Let us look at the qualitative behavior of the shell and search models which exhibits a 
critical behavior in the black hole formation. First we assume that there is some mechanism 
to prevent negative flux to emit from the shell and the shell will be static after then. It 
follows from eqs. ( |3.3|) , ( |3.11| ) and dM/dR = M/R that the orbit of the shell outside the 
event horizon (s = 1) is given by a first order ordinary differential equation 

dM ~ , „ ~ ., . k 



dR \r\ 
where / is given by eq. fl3.10| ). 



pR"(f-R" + — JP-1), (4.1) 



One immediately finds in the above equation that the value of C, which is one of the 
parameters specifying the initial condition, is irrelevant to the qualitative behavior of the 



shell, i.e. , the motion of the shell is scale invariant. This is another feature which our 
model has in common with the case of spherically symmetric scalar field collapse. The set 
T of turning points (boundary of the forbidden region) in the M-R diagram is given by the 



The line H of M = R/2 corresponds to the apparent horizon where the expansion of the 
outgoing null ray vanishes. Further, it coincides with the event horizon for the outgoing 
Vaidya spacetime. One can easily check that H is tangent to T (if u ^ 0). The shape of T 
depends on the sign of u, which is illustrated in fig. 1. If oj > then T passes the origin 
and the mass of the shell is bounded from below to be positive. If ui < the mass is not 
bounded from below. This implies that the mass of the shell can be negative at R = 0, 
which may suggest that a naked singularity occurs. We shall not concern the critical point 
between black hole and naked singularity here and assume that u is positive. The initial 
condition (Ri, Mi) can be given in the region between H and T with < R < 1, which we 
will call Region I, if the shell is not contained in a black hole from the beginning. 

One would like to consider the situation in which the shell's energy is carried away to 
infinity during the collapse. We therefore require that the mass decreases when the shell 
is collapsing. This is given by the condition p = 2r] + uj + l>0, because in Region I we 
always have \/ f 1 — 1 > / — R w > 0. Then a collapsing shell has only three possibilities: it 
will hit (1) H (R > 0), (2) T (R > 0), or (3) the origin. The case (1) implies that there is 
a formation of a black hole; the solution is supercritical. Since the local gravitational mass 
is constant along the outgoing null line, the mass Mg of the black hole, which is given by 
the limit of the Bondi mass to the future timelike infinity i + , is given by the mass of the 
shell crossing the event horizon H. The case (2) implies that the shell reaches the turning 
point and expands. But eq. ( |4.1[ ) implies that negative flux emerges when it expands so 
that the shell will stay at the turning point. The solution is subcritical. If both (1) and (2) 



points at which R = 0. From eq. ( |3.10| ) this implies / = 1 so that the mass of the shell at 



a point on T is expressed as 




(4.2) 



S 



are possible there may be a solution of the case (3), which should give the critical point. 

Whether each case really occurs can be examined by using the gradient dM/dR given 
by ( |4. 1| ) . On H one has dM/dR = 0, which implies that black hole solutions always exist. 
The gradient is always positive in Region I and negative outside Region I. Therefore, the 
shell loses/gains mass before/after it cross the event horizon. In Region I the gradient 
monotonically increases when the point is going down to T along the constant R line, 
because the partial derivative of the right hand side of (|4.1j) by M is always negative in 
Region I. For each point of H, one can check whether there is a solution with decreasing R 
which hits it by comparing the gradient of the solution and that of the line T. The difference 
is given by 

(£L-t--*('-(' + -D*) 

The point at which this is positive (negative) is the turning point from contraction to 
expansion (expansion to contraction). One finds that it changes sign once from positive to 
negative at R = Rs = (f]/(v + 1/4)) 1 /" in < R < 1 if rj > 0, and it is always negative 
otherwise. This implies that solutions of the above case (2) certainly exists if and only if 77 
is positive. Combining these conditions, we find that the critical point exists if and only if 
the both of rj and uj are positive. A few remarks are now in order. First, the condition rj > 
implies that the intrinsic force of the shell is not tension but pressure. This is a reasonable 
condition for the collapsing shell to expand. Second, all solutions with decreasing R has its 
"starting point" on T where R > Rs- The "terminal point" is either on H or on T with 
R < Rs- Third, the point R = Rs on T is the solution of the static shell. 

The critical solution M = Mq ( |4.1| ) is the one which passes the origin. It is easily found 
that Mc oc R 1+UJ + o(R 1+UJ ), in particular, the critical solution is tangent to T at the origin. 
As already said, the black hole mass Mb is given by the value of M at the intersecting point 
of the solution and M = R/2. 

The Penrose diagrams of the future half of each case are shown in fig. 2. In the critical so- 
lution the asymptotic behavior Mq oc R 1+uj implies that the invariant R a i )C dR abcd = 48M 2 /r 6 
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diverges when the shell shrink tor = 0if0<u;<2, while it is finite if lu > 2. The former 
case can be interpreted that the shell collapsed to a naked singularity in the center, while 
the latter that the shell evaporated. In both cases the curvature is always continuous in 
r > 0. 



V. APPROXIMATE SOLUTIONS 

In this section we consider the motion of the shell in some datail and obtain approximate 
solutions by an intuitive argument. We shall calculate the critical exponent by using them. 
We shall give a more rigorous argument in the next section. 

We change variables for easier manipulation of the differential equation (|4.1| ) by 



x = R" (5.1) 
y = 2MR- 1 . (5.2) 



Eq. ( fOp for R < is given by 



| = s(&>- 1 )|-^-"((; + *)'- 4 ) I- < 5 - 3 > 

The behavior of the solutions is illustrated in fig. 3. The set T of turning points is given by 
the parabolic curve y = x(2 — x). The critical solution y = yc in y-x diagram hits the origin, 
because Mc oc R l+UJ + o(R l+lJJ ) implies that yc/x — > as x — > 0. On the other hand, the 
supercritical solutions diverge at R — 0. The event horizon H is given by the straight line 
y — 1. The mass of the black hole Mb is related to the value xb of x at the point intersecting 
H as Mb oc x%. A convenient choice for initial conditions of (R, M) = (Rj,Mj) is those 
with fixed R, which correspond to initial conditions of (x,y) = (xj,yi) with fixed xj. By 
definitions of x and y one has a relation AMj oc Ayj for the initial values, where A denotes 

3' 

the difference from the critical solution. Then a power-law relation xb oc Ayj implies 
Mb oc AMj , where j3 = 13' /u. We look for the value j3' in the following. 
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Let us devide Region I into two, x > x cut and x < x cut , with x cut being a small positive 
constant. We will call the first region the linear (deviation) region. In the linear region, a 
small deviation from the critical solution is always linear, i. e. 

Ay cut cc A Vl . (5.4) 

It is because in x > x cut the rate of change of the gradient, (d/dy)(dy/dx), is bounded so 
that the deviation vector is always finite. It follows that nonlinearity of the deviation of the 
Bondi mass from the critical solution emerges only near x = 0. 

The region x < x cnt will be called the scaling region. If we ignore x compared to unity, 
eq. ( |5.3| ) reads 

d l= l -i(v-^-v{& -A\ I. (5.5) 



dx uj \ x \\x 

Note that y/x is not necessarily large compared to unity for near-critical solutions. One 
finds a scaling law in the sense that the transformation (x, y) i— > (ax, ay) does not change 
the differential equation (|5.5|), which says that if y = yi(x) is a solution, y = a^yi^ax) gives 



general solutions. In fact, one integrates eq. ( |5.5| ) to find 

fc = x( v // 2 -l + /)^ Py/f 2 -l-qf 



—J 7 

p4 —q* 



(5.6) 



where q = 2i], k is the constant of integration, and 2f = y/x in the present approximation. 
The critical solution is the one passing the origin hence the one with f = y/2x being finite 
when x —>■ 0. It is the k = solution, namely 

2 P X v\ 

VC = , I g » ( 5 - 7 ) 

Positive k corresponds to all of the other solutions. The supercritical/subcritical solutions 
are the ones with the term between the vertical lines in ( p.6|) positive/negative, because / 
should be greater/less than the critical solution f c . For small k, by setting y — 1 one finds 

k~xz^~ q . (5.8) 
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Setting x = x cut one has 

k~ (Ay cnt )tt, (5.9) 



ujq 



where one sees that the power is determined by the power of the term between the vertical 
lines in ( |5.6| ). From eqs. (p.4[), (|5.S|), ( |5.9| ) and (3 = P'/co one obtains the relation Mb oc Aff 
with the critical exponent 

(3 = — — = -. — -. (5.10) 



VI. EXACTNESS OF THE CRITICAL EXPONENT 

In the previous section we found the value of the critical exponent in a very simple 
approximation. However, it can be shown that the value actually is exact by a more rigorous 
treatment of eq. ( |4.1| ). The idea is that the cutoff scale a; cut now is dependent on the initial 
deviation of yj such that a; cut tends to zero when Ayj — > 0. Several points with different 
velocities of convergence divide the solution curve into a corresponding number of segments. 
We shall find a nonlinear relation of the deviations on each interval which eventually gives 
the value of (3. 

We rewrite eq. ( |5.3|) in terms of / = (l/2)(y/x + x) instead of y to have 

! = -^/-n/7^I-(,H>). (6.D 

Let us briefly show the existence and uniqueness of the critical solution, a solution which has 
a finite limit at x — > 0. We consider solutions of (|6.1|) in U — (0, 5) x (f CQ — (2q + l)5/uj, fco + 
(2q + 1)8 /ui) where fco = p/Vp 2 — Q 2 an d o is a positive constant (fig. 4). It can be easily 
shown that there exists 5 such that df/dx < constant < on line h : f — fco G U and 
df ' jdx — (—(2q + \)x/u) > constant > on Z 2 : / = fco — (2q + 1)x/uj G U. A solution 
cannot pass both of l\ and I2 in U. Moreover, sets of points of solutions passing l\ and of 
those of solutions passing l 2 form open sets in U. This implies that there must be a point 
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which does not cross neither li and 1% in U fl {(x, /); /co — (2g + 1)x/uj < f < fco}- So there 
must exist a solution f = fc which satisfies 

(2g + l)x 

/co < to < /co (6.2) 

in < x < 5 and has a limit /<? — > fco as x ^ 0. 

Let us define A by A = f — fc- We consider solutions with nonnegative A. The equation 
for A is given by 

' /1UA --(pF(fJc)-q), (6.3) 



din a; uo 



where 



F(g, b) = F(b, a) = , Q + (6.4) 

for a, b > 1. 

We show some inequalities derived from eq. ( |6.3j ). F(a, b) is a decreasing function of a 
and b so that 

F(a,&)>F(c,c) = ^==>l (6.5) 

V(T — 1 

holds for c> a,b. Since / > /c? > 1, if / < c on interval [x min , x max ] then 

din A 1 ( pc \ p-q 

dlnT " ~u ~ 7 ^~ (6 ' 6) 

holds on the interval. Its integration on interval [a: mm , x max ] gives 

A-(x m in) ^ / Smax \ ~ ( Vc 2 -! ~ 9 ) (6 7) 



•aim y . . / ■"max \ - V \ - l 

AfXrnnv) V X 



mm / 

£ 1 

max 



> (^) ~ . (6.8) 

In particular, A increases as x decreases and diverges as x — > 0. This implies that the critical 
solution is unique. 

By this and a similar argument for negative A, one finds that all solutions / > fc go to 
oo as x — > 0, while ones f < fc to — oo (though, in fact, / < 1 is forbidden). This tells that 
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the point (0, fco) is a separatrix. One expects that the deviation from the critical solution 
becomes nonlinear near the separatrix. 

If d < fc < f < c on the interval [x m i n , x max ] with c — d being sufficiently small then 

F(/ , /o) <_J rT <_^ T + _f^ (6 . 9) 



holds, where the second inequality is because c/ a/c 2 — 1 is a convex function of c. This gives 
cHnA 1 / »c p(c — d) \ 

dl^ £ "(7fe- 9+ (^T)37ij' <610) 
and its integration yields 

A(x min ) < / x max \ £ (7fe~ g + (J^)!^ ) (6 11) 



If fc < < c < / on [x min , x max ] with sufficiently large c and fixed d it is easily shown that 

F{f,f c )<l + — (6-12) 
c 

holds on [xmm, x max ], so that one has 

A(x ■ ) fx \b{p-i+ s ^) 

< . (6.13) 



-Trim 



A(3? max ) V X n 

For notational simplicity we adopt the convention that subscript i(= 1,2,3,...) denote 
the value at x = X{ and subscript denotes the value at x = 0. Let x± = xj and A x = e (fig. 
5). Let us define x 2 = | lne| _1 . From eq. ( |6.8|) one has 



p-i 



At ~ \Jlne, 
From eq. (|6.3|) one has 

din A 1 / p{2f c + A) 



dlnx u VvT^T + v^/c^T V 

>-^f-^=-g + -^—|, (6.15) 
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PA(l + QA)d\nx 
where P and Q are positive constants. Integration on [x 2 ,£i] gives 



dA 

> -1, (6.16) 



1+ ^>m p , (6.i7) 



1 + ofc 1 



hence 



In e| _1 \ „ / , /I lneT 1 ^ 
- Qe 1 



< 



X\ ) \ \ X\ 

lnel" 1 \ P 1 / 1 In e| _1 \ 



X\ j 2\ xi 

\ p 
X\ \ 



-i 



(6.18) 



y | In 

where the inequality between the second and the third lines holds for sufficiently small e. 
Combining inequalities ( p. 14 ) and ( |6.18| ) one has 



A 2 ~ Ai ~ e. (6.19) 

Here A ~ B means that A and B are the same in power of e up to logarithm. More 
precisely, A ~ B holds if there exist positive constants, i.e. numbers independent of e, P, P', 
and constants Q, Q' such that 

P|lne| Q < - < P'\\ne\ Q ' (6.20) 

holds. 

Let us define £3 such that A 3 = | lne| _1 . From A 2 ~ e one has £3 < x 2 for sufficiently 
small e. Inequality (|6.2|) and the definition of 23 implies that 



, 2( 1+ !, 1 1-1 ^ f ^ f 
fco 1 lne| < f c < fco, 



fc<f<fco + line)- 1 . (6.21) 
It follows from eqs. ( |6.21|) and (|6.7p that, if e is sufficiently small, 
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~2 ~ VA3- 



(6.22) 



holds for some constant P. From inequalities ( |6.21| ) and ( |6.11| ) 

/A\ ^T+P'llnel" 1 

x 2 ~ VA 3 , 



(6.23) 



holds for some constant P' . It follows from these inequalities and the orders in e of 22, A 2 
and A3 that 



ujq 
2 _ 2 



(6.24) 



because the term P'\ lne| 1 in the exponent causes only logarithmic corrections in ( |6.23| ). 
Let us define X4 such that A4 = | lne|. Monotonicity of A implies 



x 4 < x 3 . 



(6.25) 



Eq. (|6.8|) implies 



From these inequalities one has 



%i > /Aj^ p- 
x 3 ~ VA 4 , 



Lxjq 

x 4 ~ x 3 ~ 



(6.26) 



(6.27) 



Let us define x 5 = xb- The definition of xb implies 

1 / 1 



It can be easily checked that 



A 5 = h x 5 - f CB . 

2 \x 5 



(6.28) 



fc < fco < A 4 = |lne| < /. 



(6.29) 



From eqs. ( |6.29 ) and (O) one has 



A 5 /x 4 
A 4 ~~ Vx 5 



(6.30) 
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This implies 



x 5 > Px 4 p ~ q A A w (6.31) 

for some positive constant P, because A 5 < 1/2x5 f° r sufficiently small e and p — q — ou + 1. 
One has from ( |6 .13 ) that 



A 5 /x 4 \z;{p-i + W7\) 

< — V ' (6.32) 



A 4 \x 5 
for a constant P' . These imply that 

x 5 ~ x/- q ~ eF%. (6.33) 

As a result, we have 

x B ~Ai^, (6.34) 
which proves that the critical exponent (3 is equal to qj (p + q). 



VII. SUMMARY 

We constructed an analytic model which exibits a critical behavior of black hole forma- 
tion, namely, collapse of a thin shell coupled with outgoing null fluid. In that model the 
dynamics of the shell is described by an ordinary differential equation. We investigated the 
dynamics of the shell and the evolution of the Bondi-Sachs mass. It was seen that the devi- 
ation of the Bondi-Sachs mass from that of the critical solution becomes nonlinear only near 
the M = and near the event horizon. With suitable choice of dymamical variables / and 
x, this is, as expected in critical phenomena, understood as follows: the point corresponding 
to the final state of the critical solution is represented as a separatrix, and the nonlinear 
relation AMg oc AM/' 3 for the supercritical solutions occurs near that point. We calculated 
the value of the critical exponent (3 analytically to find (3 = 2/7/(477 + u + 1), where 77 is 
a parameter in the equation of state of the thin shell and uj specifies the strength of the 
coupling of the thin shell and the outgoing flux. 
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VIII. DISCUSSION 



It is remarkable that in our model the critical exponent (3 is dependent on the constants 
T) and u specifying the system. The exponent (3 can take the values in the range < (3 < 1/2 
for i] > and uo > 0. This suggests that the critical exponent is not so universal in the 
extended sense in Sec. [I| There may be many different universality classes, and the critical 
exponent (3 may serve as an index specifying them. 

There will be a question why the numerically calculated cases seem to fall into the same 
universality class with (3 ~ 0.38, even if we admit the existence of many universal classes. 
One may think that the massless property of the matter is crucial. In our case, even if 
we require that the shell consists of radiation fluid having a trace free equation of state, 
r\ — 1/2, then (3 is equal to 1/(3 + uj) and cannot exceed 1/3. This may suggest that the 
massless property is not so important. Another possibility is that self-similarity is crucial. A 
support for this possibility is that the critical solution in our model does not have spacetime 
self-similarity. If the near-critical solutions necessarily fall into a self-similar spacetime in 
the systems with smoothly distributed matter, then types of self-similar spacetimes may 
determine the universality classes. 

The universality in the initial value is not verified in our model because the solutions 
are parametrized by a single parameter Mj (or Rj), apart from an irrelevant parameter C . 
In our model, as said above, the nonlinearity of the deviation of the mass from the critical 
solution emerges only near the event horizon. If the situations are similar in the numerically 
investigated collapses, the universality in initial data can be naturally understood. Namely, 
almost all modes or fluctuations in initial data damp away so that the spacetime fall into 
a unique one by the time that the event horizon forms, and the critical exponent depends 
only on that final portion of the spacetime. This may also explain the universality in the 
extended sense if the spacetimes fall into the same one in the last stage of the collapse in 
several different systems. 

To understand the universality it may be instructive to parametrize the solution space 



18 



differently, for example, pi = CMi and P2 = C~ x Mi. Then the critical exponent is, of 
course, the same for the both parameters. In this point of view, one can say that the 
universality on initial data is verified in our model and that the choice of the parameters 
Mj and C was good one in the sense that it separates the relevant and the irrelevant mode 
to the criticality. In more realistic collapses, where there are many degrees of freedom, it 
is expected that there is one (or a few) relevant mode corresponding to Mj and are lots of 
irrelevant modes corresponding to C. 

We remark that our model has a peculiarity that if we do not assume some mechanism 
to protect the energy condition and allow negative flux, the subcritical solutions will have 
a reexpansion of the shell and a formation of a black holes of finite mass [§. Though the 
above assumption seems somewhat artificial, we believe that in more natural models the 
mechanism of the emergence of nonlinearity is similar to that of the model presented here. 

Some authors have argued black hole evaporation by imitating the spacetimes with the 
Vaidya spacetime or one with classical scalar field (e.g., |10| , |TT|| ). If these point of view have 
something in common with the real machanism, some critical phenomena may be found in 
black hole evaporation and would be the key to understand such spacetimes. 
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FIGURE CAPTIONS 



Fig. 1. The evolutions of the Bondi — Sachs mass and the radius of the shell are illus- 
trated. The line H corresponds to the apparent horizon and the line T is the set of states 
with R = 0. Region I is between H and T, where we can set initial data. The thick line 
C with an arrow denotes the critical solution which hits the origin. The thin line with an 
arrow above C is a supercritical solution which cross H. The value M B of the mass function 
is the mass of the formed black hole. The thin line with an arrow below C is a subcritical 
solution. 

Fig. 2. The Penrose diagram of (a)the supercritical solution, (b)the critical solution, 
and (c)the subcritical solution, where the double lines denote singularities. In (b) the r = 
line above the point of contraction of the shell should be replaced by double line. 

Fig. 3. The evolution in variables x and y. The line H corresponds to the apparent 
horizon and the line T is the set of states with R = 0. Region I is between H and T, where 
we can set initial data. The thick line C with an arrow denotes the critical solution which 
hits the origin. The thin line with an arrow above C is a supercritical solution which crosses 
H . The coordinate of the point intersecting H gives the mass of the formed black hole. 
The thin line with an arrow below C is a subcritical solution. 

Fig. 4. The behavior of the solutions in variables x and /. The curves are the solutions. 
The arrows on li and li indicate gradients of the solutions at those points. There must exist 
some solution f — fc which has a finite limit lim^o / = /co- 
Fig. 5. A figure which illustrates the procedure of finding the exact value of the critical 
exponent f3. 
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